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ABSTRACT 

We present preliminary results from a programme to obtain multi-object near- 
infrared spectroscopy of galaxies at redshifts 0.7 < z < 1.5. We are using the in- 
strument CIRPASS (the Cambridge Infra- Red PAnoramic Survey Spectrograph), in 
multi-object mode, to survey Ha in galaxies at z ~ 1. We aim to address the true star 
formation history of the universe at this epoch: potentially the peak period of star 
formation activity. Ha is the same star formation measure used at low redshift, and 
hence we can trace star formation without the systematic uncertainties of using differ- 
ent calibrators in different redshift bins, or the extreme dust extinction in the rest-UV. 
CIRPASS has been successfully demonstrated in multi-object mode on the AAT and 
WHT. Here we present preliminary results from one of our fields, the Hubble Deep 
Field North, observed with the WHT. With 150 fibres deployed over an unvignetted 
field of ~ 15arcmin, we have several detections of Ha from star forming galaxies at 
0.8 < z < 1.0 and present spectra of the seven brightest of these. By pre-selecting 
galaxies with redshifts such that Ha will appear between the OH sky lines, we can 
detect star formation rates of 5 /if 2 Mq yr" 1 (5 a in 3-hours, Qm = 0.3, = 0.7 ). It 
appears that star formation rates inferred from Ha are, on average, a factor of more 
than two higher than those based on the UV continuum alone. 
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1 INTRODUCTION 

A central problem in observational cosmology is determin- 
ing at which epoch the majority of stars formed. This has 
important implications for models of galaxy formation and 
evolution. Despite recent progress in this area, the star for- 
mation history of the universe is still a topic of intense de- 
bate. There is substantial evidence that the star formation 
rate was much higher in the recent past, compared with 
the current epoch, rising steeply to z ~ 1 (e.g. Lilly et al. 
1996, Tresse et al. 2002, Hippelein et al. 2003). However, it 
is still unclear whether at higher redshifts the star formation 
density declines, plateaus or perhaps continues to slowly in- 
crease. Most quantitative attempts to measure the global 
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star formation history have suffered from having to use dif- 
ferent indicators of star formation in various redshift bins 
(Madau et al., 1996), redshifted into the optical. These var- 
ious indicators not only have uncertain relative calibration 
but are also affected differently by dust extinction. To make 
a reliable comparison, the same tracer of star formation must 
be used at high redshift as that used locally. The Ha emis- 
sion line is a good tracer of the instantaneous star formation 
rate, and has been widely used in surveys at low redshift. It 
is particularly suitable as it is relatively immune to metallic- 
ity effects and is much less susceptible to extinction by dust 
than the rest-UV continuum and Lyman-a (which is also 
selectively quenched through resonant scattering) . It is also 
important to use Ha to calibrate secondary indicators of star 
formation such as [OII]3727A, which are often used at red- 
shift one. However, tracing Ha to early epochs forces a move 
to the near-IR at z > 0.6. The recent advent of good infra- 



2 M.Doherty et al. 



red spectrographs on large telescopes has facilitated rapid 
advancement in this area, with observations of a few tens of 
galaxies at z ~ 1 (e.g. Tresse et al., 2002, Glazebrook et al., 
1999) and z ~ 2 (Erb et al., 2003). However, until recently, 
near-infrared spectroscopy has been restricted to long-slit 
work and building samples using single object spectroscopy 
is inefficient in terms of telescope time. The small statis- 
tical samples obtained result in large uncertainties in the 
global properties of galaxies at z ~ 1. There has been suc- 
cessful "multi-object" spectroscopy through slitless surveys 
in the near-infrared from space (the HST/NICMOS survey 
of Yan et al. 1999), but such an approach has poor sensi- 
tivity because of the high background. It is only now that 
true multi-object infrared spectroscopy is possible from the 
ground, using either a slitmask approach (e.g. IRIS 2 and 
FLAMINGOS) or a fibre-fed spectrograph. In this paper 
we present the first successful demonstration of multi-object 
near-infrared spectroscopy of high-redshift galaxies. 

We have used our fibre-fed CIRPASS spectrograph 
(Cambridge InfraRed PAnoramic Survey Spectrograph, 
Parry et al. 2000) to measure the star formation rates of 
a sample of galaxies in Hubble Deep Field North (HDF-N) 
(Williams et al., 1996). This is the initial stage of a larger 
survey to address the true star formation history of the uni- 
verse at redshifts z = 0.7 — 1.5, through Ha measurements 
of several hundred galaxies. CIRPASS can operate with a 
150 fibre multi-object bundle, with the ability to simultane- 
ously observe 75 object+sky targets. It thus offers a huge 
multiplex advantage, compared with surveys utilising single 
object spectroscopy - of which the largest to date is the work 
of Tresse et al. (2002), which surveys a total of 33 galax- 
ies. Controlled selection of our targets ensures an accurate 
knowledge of the completeness of our survey. The selection 
function will be discussed in a future paper (Doherty et al., 
2004, in prep.). 

In this paper, we focus on the Ha detections from our pi- 
lot study of galaxies at z ~ 1 in the HDF-N with CIRPASS. 
The layout of the paper is as follows: in section 2 we describe 
the instrument and our observations and in section 3 we de- 
tail the data reduction techniques employed. In section 4 we 
present results for a sub-sample of objects, which constitute 
our brightest detections of Ha (with S/N > 5 a), and use 
the results to derive the star formation rates in these red- 
shift one galaxies, comparing these with that inferred from 
the rest-frame UV (2400A) continuum. Our conclusions are 
presented in section 5. 

In this paper we adopt the standard "concordance" 
cosmology of Qm = 0.3, £Ia = 0.7, and use /170 = 
// /70kms" 1 Mpc" 1 . AB magnitudes (Oke & Gunn, 1983) 
are used throughout. 



2 OBSERVATIONS 

We performed multi-object spectroscopy of galaxies in the 
vicinity of the HDF-N with CIRPASS, a near-infrared fibre- 
fed spectrograph operating between 0.9 and 1.67 /jm. The 
upper cut-off is set by a blocking filter which reduces the 
thermal background. These HDF-N observations were un- 
dertaken at the cassegrain focus on the 4.2 m William Her- 
schel Telescope (WHT) in La Palma. The instrument was 
used in multi-object mode, with 150 fibres of 1.1" diam- 
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Figure 1. Spectra for the 7 brightest Ha detections in our sam- 
ple. The expected position of Ha at the optical spectroscopic red- 
shift is marked with a dashed line. The sky spectrum is overlaid 
in dotted lines. The emission lines all fall well between sky lines. 
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eter on the sky (comparable to the expected seeing con- 
volved with typical galaxy profiles) . These fibres are deploy- 
able over an unvignetted field of diameter 15'. We used the 
FOCAP mechanism developed by the AAO as an interface 
to the telescope, with holes drilled in a brass plate to hold 
fibres in position. We acquired the field using six guide bun- 
dles, each containing seven closely packed fibres, centred on 
bright stars. We thus centroided the plate to an accuracy 
better than 0.5", that is, less than half the fibre size. 



A Hawaii 2K detector was used, and a grating of 
831 lines/mm, producing a dispersion of 0.95A/pixel. The 
FWHM of each fibre extends over 2.7 pixels in both the 
spatial and spectral domain. The wavelength coverage was 
1726A, covering most of the J-band, and the grating was 
tilted to set a central wavelength of A c = 1.25/im. The re- 
solving power was R = X/ AXfwhm ~ 5000. At this res- 
olution the background is very dark between the OH sky 
lines in the J- and H-bands, with only about 10% of the 
array contaminated by skylines. By targeting galaxies with 
redshifts which place the Ha emission lines between the sky- 
lines we become very sensitive to line emission - indeed, we 
are limited by the instrument background (although it is 
cooled to — 42° C) rather than the sky. We targeted 65 ob- 
jects using pairs of fibres offset by 6', nodding the telescope 
between the A and B positions such that the fibre pairs al- 
ternately received light from the target object and sky. The 
2D data frames are then subtracted directly from each other 
before extracting the spectra, so that sky subtracted off has 
been observed down the same fibre as the object. Although 
the sky background varies a lot on timescales equivalent to 
our exposure time, we achieved reasonable first order sky 
subtraction. Some low throughput fibres were not allocated 
positions. A total of 12 hours integration time was obtained, 
over 3 nights, with individual integrations of 30 or 40 min- 
utes per pointing. The array was read out non-destructively 
every 10 minutes, with three or four such loops per pointing 
to allow for cosmic ray rejection. There were 10 reads per 
each 10 minute loop to reduce the readout noise, ensuring 
that we achieved background-limited observations. 



The majority of our targets in this field were selected 
from the redshift survey of the HDF-N and "flanking fields" 
carried out by Cohen et al. (2000), which is > 92% com- 
plete to Rab — 24 for objects in the HDF-N proper and 
Rab = 23 in the flanking fields. We selected all Rab < 24.5 
targets within the redshift range 0.73 < z < 1.0 which fell 
within our 15' unvignetted field of view. A handful of ob- 
jects were removed because of the limit on the fibre spacing 
(a minimum of 20"). We redefined the coordinates of our 
targets based on the GOODSvl.O system (Giavalisco et al., 
2003), for consistency in our astrometry between the targets 
and the alignment stars. The completeness of our sample is 
discussed in Doherty et al. (2004, in preparation). In this pa- 
per we present observations of the seven sources in HDF-N 
where Ha was detected at more than 5a. These sources are 
all in the Cohen et al. (2000) flanking field, with Rab < 23.0. 



3 DATA REDUCTION 

Data reduction was performed under IRAF with the CIR- 
PASS package 1 . For each multiple read exposure, the aver- 
age of each loop was subtracted from the average of the next 
loop of non-destructive reads to form sub-integrations of 10 
minutes each at the same pointing. These sub-integrations 
were then averaged together using the crreject algorithm 
to reject cosmic ray strikes (using limits of 3<r according to 
gain and readnoise). The frames were then beam-switched 
by taking A-B pairs (providing first-order sky subtraction). 
We then subtracted of the average residual bias, determined 
from the unilluminated region of the array (the bias level 
of infrared arrays is known to float with time). The spec- 
trum in each fibre was then optimally extracted according 
to the prescription of Johnson et al. (2002), that is, the fibres 
were profile fitted in triplets (central fibre and two adjacent 
fibres) and the pixel values weighted and summed. This re- 
moves the contribution from the two adjacent fibres to the 
central fibre. This process produced positive (object-sky) 
and negative (sky-object) fibre pairs. 

Dome flat fields were also optimally extracted and used 
to correct variations in the fibre throughput and the spectral 
response of pixels on the array. Wavelength calibration was 
performed using an argon arc, fitting 30 lines with a third 
order polynomial, resulting in an rms dispersion of O.lA. The 
spectra were then rectified to lie on a common wavelength 
scale. 

In order to remove sky residuals which did not fully 
subtract out in the beam-switch, a low order polynomial 
was fitted and subtracted along each rectified column of the 
flat-fielded, extracted spectra. This removed residuals due 
to temporal variation of the night sky spectrum. Flux cal- 
ibration was carried out using observations of a J = 10.68 
magnitude star from 2MASS which was placed on the same 
HDF-N plate as our targets. As it was thus observed si- 
multaneously, this also corrects for temporal changes in the 
seeing, which fluctuated a lot over the course of the obser- 
vations. As a consistency check we looked at 6 bright stars 
placed on the plate in another field from the survey. The 
fluxes were consistent to within 3%, comparing the counts in 
our CIRPASS spectra with those predicted from the 2MASS 
J— band magnitudes. 

In an average fibre, between skylines, for a spectrally 
unresolved line (< 60km s _1 ) we achieve a sensitivity of 7.2 x 
10~ 17 crg s _1 cm" 2 at 5a in 3 hours. However the emission 
lines we detect are typically broader than this (~ 100 — 
250 km s _1 , Table 1). 



4 RESULTS 

In this paper we analyse a sub-sample of our HDF-N data: 
those galaxies with strongly detected (greater than 5a) Ha 
emission. These emission line spectra are shown in Figure 1. 
The spectrum of the night sky is overlaid in dotted lines, 
giving a clear indication of the positions of skyline residuals, 
where the noise is significantly higher due to the enhanced 
background counts. These residuals do not adversely affect 
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the data as the emission lines shown here all fall well between 
skylines. The dashed line shows the expected position for 
Ha, given the optical spectroscopic redshift. The centroid 
of Ha is accurate to 2 pixels (~ 2A), so the error on the 
redshift is less than 0.001. There is agreement at this level 
with the optical redshifts given in Cohen et al. (2000) and 
which are listed in Column 5 of Table 1. 

We accurately measured the integrated line fluxes by 
measuring between zero power points (Table 1), having first 
subtracted off a continuum fit derived from regions adja- 
cent to the emission line and unaffected by skylines. We 
checked our flux measurements by also fitting gaussian pro- 
files to the lines, finding consistent values. The gaussian fits 
also gave the full width at half maximum (FWHM) for each 
line (Table 1). After subtracting the instrumental resolu- 
tion in quadrature from the galaxy line widths, we find 
velocity FWHMs for these galaxies in the range ~ 100 — 
250 km s -1 , equivalent to am — 40 — 100 km s^ 1 . We note 
that some of this velocity width may be due to galactic rota- 
tion in the extended galaxies (galaxies J1236175+6214027, 
J1237084+6215150 and J1237166+6210424 are the most 
spatially extended and also exhibit the broadest lines). 

The Ha luminosity of a galaxy is directly proportional 
to the ionising flux from massive stars, which drops off 
very quickly, about 20 million years after star formation 
ceases. Ha thus traces the instantaneous star formation 
rate, whereas the UV luminosity evolves in time with the 
changes in stellar population and continues to rise even af- 
ter star formation has ended. Since the Ha flux is propor- 
tional to the number of OB stars, in order to extrapolate 
to a total star formation rate (SFR) it is necessary to as- 
sume an initial mass function (IMF), and the conversion 
from Ha flux to SFR is quite sensitive to the IMF assumed 
(see discussion in Glazebrook et al. 1999). We use Kenni- 
cutt's (1998) conversion which assumes a Salpeter IMF with 
0.1 Mq <M < 100 Mq : 

SFR(M Q yr- 1 ) = 7.9 x 10~ 42 L Ha (erg s" 1 ) (1) 

The derived SFRs of the galaxies are shown in Table 2, 
and in Figure 3 we show B, V, i' and z' band images from 
HST/ACS taken from GOODS vl.O (Giavalisco et al., 2003). 
We performed photometry on these images with 1" aper- 
tures, for consistency with our infra-red fibre size. In most 
cases the 1.1" fibres enclose most of the B-band light. How- 
ever, for three of the sources (galaxies J1236175+6214027, 
J1237084+6215150 and J1237166+6210424) some fraction 
(< 20%) is missed. We reiterate that we have compensated 
for the effects of seeing by placing one of the fibres on a 
bright star in the field (thus observing the star simultane- 
ously with the galaxies). This allows us to correct the flux 
calibration for seeing dependant aperture losses. 

We used the B band magnitudes (4500A) to calculate 
rest frame UV (2400A) flux densities and corresponding star 
formation rates. These are also shown in Table 2. For consis- 
tency, we also use the conversion given in Kennicutt (1998): 

SFR(M yr" 1 ) = 1.4 x 10~ 28 L„(erg s" 1 Hz" 1 ) (2) 

assuming the same IMF as used when deriving SFRs from 
the Ha flux (Equation 1). The UV SFR relation also assumes 
continuous star formation over ~ 10 8 years. 

Figure 2 shows the SFRs for each galaxy calculated us- 
ing the UV luminosity density and the Ha flux: those cal- 



culated from the UV luminosity densities are a factor of 
two lower on average, than those from Ha. This is proba- 
bly due to the differential effect of dust extinction in the 
redshift one galaxies between A rcst « 2400A and 6563A. 
This is consistent with results obtained by Glazebrook et 
al. (1999), Tresse et al. (2002), and Yan et al. (1999) who 
all find SFR(Ha)/SFR(UV) ratios of around 2-3. 

There is no uniform correction that can be applied to 
the SFR estimated from UV flux in order to derive the 
true star formation rate: the amount of extinction varies 
in each object due to inherent differences in physical prop- 
erties of the galaxies (e.g. (Sullivan et al., 2001)). As can be 
seen in Figure 2, the statistical errors in our Ha flux mea- 
surements are not great enough to account for the scatter, 
which instead is attributable to different dust extinctions in 
our galaxy sample. Ha provides a more robust indicator for 
SFR, giving an independent estimate that is much less af- 
fected by dust obscuration. The ratio of the total SFRhc* to 
SFRyv for this sub-sample is 2. This fraction is potentially 
biased high, as there is a selection effect in this sub-sample 
due to the fact that we are only including the sources with 
strong detections of Ha: we will address the full i?-band 
magnitude-limited sample in a future paper. 

We now estimate the true star formation rate for 
each galaxy, assuming that the Kennicutt (1998) relations 
for Ha and UV SFRs should yield the same extrinsic 
SFR in the absence of extinction (i.e. assuming a Salpeter 
IMF and continuous star formation). Furthermore, we as- 
sume a Calzetti dust extinction law, appropriate for star- 
burst galaxies (Calzetti, 1997). For each galaxy we calcu- 
late the reddening value E(B-V) which gives the observed 
SFR(Ha) /SFR(UV) ratio. In Table 2 we list the reddening 
values derived and use these to compute the total dust cor- 
rected star formation rates. These are typically 50% higher 
than those based purely on Ha, corresponding to an aver- 
age E(B — V) = 0.16. Hence surveys based solely on UV 
(e.g. 2400A) continua (without dust correction) will under- 
estimate SFRs at redshift one by factors of « 2 — 11. 



5 CONCLUSION 

We have performed multi-object, near-infrared spectroscopy 
of z ~ 1 galaxies in the Hubble Deep Field North, using 
CIRPASS-MOS. These observations are part of an on-going 
survey to trace star formation rates at redshift one. We have 
presented our brightest detections of Ha from three hours 
of observing time, as a demonstration of the success of this 
technique. We have shown that we can detect Ha at suffi- 
cient signal to noise to obtain a good estimate of star forma- 
tion rates in these galaxies. By pre-selecting galaxies with 
Ha redshifted between the OH sky lines, we can detect star 
formation rates of 5 /if 2 Mq yr" 1 (5 a in 3- hours). The star 
formation rates obtained are higher than those estimated 
from UV continuum by a factor of about two, due to dust 
obscuration in the UV. This is consistent with previous work 
in this area. We have another ~ 60 sources in this field, at 
known redshifts, therefore by stacking the spectra we will be 
able to obtain a global star formation rate for z ~ 1 galaxies 
in the vicinity of the HDF-N. Between this and other fields 
we hope to build up a sample of several hundred galaxies 
which will allow us to determine the Ha luminosity function 
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id 


name 


RA 


Dec 




Zopt 




fwhm 


Vel. FWHM 


flux xl0~ 16 






(J2000) 


(J2000) 


(A) 




(A) 


km s _1 


erg s" 1 cm~ 2 


(a) 


J1236175+6214027 


12 36 17.536 


+62 14 02.70 


11924 


0.818 


0.8169 


* 


* 


3.60 + 5.58 


(b) 


J1237063+6215185 


12 37 06.293 


+62 15 18.50 


12078 


0.84 


0.8404 


5.9 


127 


3.31 ± 0.38 


(c) 


J1237084+6215150 


12 37 08.381 


+62 15 15.04 


12074 


0.839 


0.8398 


11.1 


259 


3.01 ± 0.54 


(d) 


J1237087+6211285 


12 37 08.659 


+62 11 28.52 


12520 


0.907 


0.9077 


5.8 


124 


2.54 + 0.34 


(e) 


J1237141+6210448 


12 37 14.141 


+62 10 44.78 


11952 


0.821 


0.8212 


4.8 


97 


2.27 + 0.36 


(0 


J1237166+6210424 


12 37 16.631 


+62 10 42.36 


11949 


0.821 


0.8207 


6.6 


146 


1.78 ± 0.40 


(g) 


J1237167+6213105 


12 37 16.716 


+62 13 10.54 


12465 


0.898 


0.8993 


5.5 


116 


2.63 ± 0.31 



Table 1. Measured properties of the Ha emission lines in our > 5<x detections in the HDF-N. The optical redshifts are from Cohen et 
al. (2000). The Ha line in J1236175+6214027 is not well fit by a gaussian profile, so no FWHM is reported. 



id 


Rab 


SFRhc 


Bab 


flux density L v 


SFRyv 


ratio 


E(B-V) 


SFR tot M yr" 1 




mag 


Mq yr _1 


mag 


10 28 erg s" 1 Hz- 1 


M© yr- 1 


SFRhc/SFRc/v 


mag 


(dust corrected) 


(a) 


21.73 


9.0 


23.89 


1.76 


2.5 


3.6 


0.30 


28.7 


GO 


22.22 


8.9 


22.89 


4.70 


6.6 


1.3 


0.07 


11.9 


(c) 


21.62 


8.0 


23.44 


2.81 


3.9 


2.0 


0.17 


15.4 


(d) 


22.23 


8.3 


23.23 


4.00 


5.6 


1.5 


0.10 


12.2 


(e) 


22.32 


5.7 


23.32 


3.01 


4.2 


1.4 


0.09 


8.1 


(f) 


21.59 


4.5 


23.85 


1.85 


2.6 


1.7 


0.13 


7.5 


(g) 


22.72 


8.4 


23.96 


2.01 


2.8 


3.0 


0.25 


22.1 



Table 2. Identification, Rab magnitude (Cohen et al., 2000), SFR derived from Ha, HST B-band magnitude (1" aperture) and 
corresponding flux density and SFR ul) for our > 5cr detections. 
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Figure 2. Comparison of SFRs obtained from UV continuum 
flux at 2400A versus Ha flux for the individual galaxies. The SFRs 
derived from UV luminosity arc consistently underestimated. The 
solid line has a gradient of 2.11 and represents the line of best fit 
to the data (using a least squares fit). The errors in the SFR(Ha) 
values reflect the errors in the flux calibration. We assume ~ 5% 
error in the UV flux, which is dominated by errors in positioning 
the fibres (i.e. mirrored in the exact positioning of the aperture 
for the Bab magnitude.) 

at z ~ 1 and hence address the true star formation rate at 
this important epoch. The success of our pilot multi-object 
spectroscopic survey bodes well for future instruments, such 
as FMOS on Subaru. 
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